Anderson localization and nonlinearity in one dimensional disordered photonic lattices 
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We experimentally investigate the evolution of linear and nonlinear waves in a realization of 
the Anderson model using disordered one dimensional waveguide lattices. Two types of localized 
eigenmodes, flat-phased and staggered, are directly measured. Nonlinear perturbations enhances 
localization in one type, and induce derealization in the other. In a complementary approach, we 
study the evolution on short time scales of 5- like wavepackets in the presence of disorder. A transition 
from ballistic wavepacket expansion to exponential (Anderson) localization is observed. We find an 
intermediate regime in which the ballistic and localized components coexist while diffusive dynamics 
is absent. Evidence is found for a faster transition into localization under nonlinear conditions. 



The propagation of waves in periodic and disordered 
structures are at the foundations of modern condensed- 
matter physics. Anderson localization is a key concept, 
formulated to explain the spatial confinement due to dis- 
order of quantum mechanical wavefunctions that would 
spread over the entire system in an ideal periodic lat- 
tice ll|, 0, H, 3- Although Anderson localization was 
studied experimentally, the underlying phenomena - the 
emergence of localized eigenmodes and the suppression 
of wavepacket expansion - were rarely observed directly 
@, Q. Instead, localization was usually studied indi- 
rectly by measurements of macroscopic quantities such 
as conductance Q, backscattering @, 0] and transmis- 
sion jjJGjj. 

An interesting issue concerns the effect of nonlinear- 
ity on Anderson localization. Nonlinear interactions be- 
tween the propagating waves and nonlinearly accumu- 
lated phases can significantly change interference proper- 
ties, thus fundamentally affecting localization. The the- 
oretical study of the nonlinear problem advanced using 
several approaches: the study of the transmission prob- 
lem [ll[ ; the study of the effect of nonlinear perturbations 
on localized eigenmodes [12j; and the study of the effect 
of nonlinearity on wavepacket expansion in the presence 
of disorder [13[. Only a few experiments were reported 
[5[. Recently, optical studies enabled the study of wave 
evolution in nonlinear disordered lattices [3, Qj3, Gil ? us ~ 
ing a scheme discussed in In particular, Schwartz 

et al [16[ reported the observation of Anderson local- 
ization of expanding wavepackets in 2D lattices. 

In this work we investigate directly linear and nonlin- 
ear wave evolution in one dimensional (ID) disordered 
photonic lattices, using two different approaches. In the 
first part of this work, all the localized eigenmodes of a 
weakly disordered lattice are selectively excited. Non- 
linearity is then introduced in a controlled manner, to 
examine its effect on localized eigenmodes. The second 
part of this work presents a study of the effect of dis- 
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FIG. 1: (color online), (a) Schematic view of the sample used 
in the experiments. The red arrow indicates the input beam. 
(b)-(d) Images of output light distribution, when the input 
beam covers a few lattice sites: (b) in a periodic lattice, (c) 
in a disordered lattice, when the input beam is coupled to a 
location which exhibits a high degree of expansion, and (d) 
in the same disordered lattice when the beam is coupled to a 
location in which localization is clearly observed. 



order on the evolution of J-like wavepackets (single site 
excitations). A transition from free ballistic wavepacket 
expansion to exponential localization is observed, as well 
as an intermediate regime of coexistence. We then mea- 
sure the effect of nonlinearity on this process. 

Our experimental setup is a one-dimensional lattice 
of coupled optical waveguides patterned on an AlGaAs 
substrate pj| [2(j, illustrated in Fig. la. Light is in- 
jected into one or a few waveguides at the input and can 
coherently tunnel between neighboring waveguides as it 
propagates along the z axis. Light distribution is then 
measured at the output (see for example Fig. l(b)-(d)). 

The equations describing light dynamics in these struc- 
tures are identical (in the linear limit) to the equations 
describing the time evolution of a single electron in a lat- 
tice under the tight binding approximation [19], i.e. a set 
of coupled discrete Schrodinger equations: 



.dU n 



f3nU n +C n , n±1 (U n+1 + f/„_ 1 )+7|C/ n | 2 C/„ (1) 



. Here n = 1, N where ./V is the number of lattice sites 
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FIG. 2: (color online). Measurements of localized eigenmodes of a disordered lattice, (a) Calculated eigenmodes and eigenvalues 
of the weakly disordered lattice used in the experiments. The band of eigenvalues deviates slightly from the cosine-shaped band 
of a periodic lattice. Localized modes are formed, associated with eigenvalues near the edges of the band (insets 1,2), while 
modes near the band center remain extended (inset 3). (b) Measurements of pure flat-phase Anderson localized modes. Panels 
show a comparison between measurements (blue) and the corresponding calculated eigenmodes of the lattice (red), (c) Same 
for staggered localized eigenmodes. In all cases no fitting procedures are used. 



(waveguides), U n is the wave amplitude at site n, (3 n 
is the eigenvalue (propagation constant) associated with 
the n'th site , C n?n ±i are the tunnelling rates between 
two adjacent sites, and z is the longitudinal space coor- 
dinate. The last term in Eq. (1) describes the nonlinear 
dependence of the refractive index on the light intensity, 
where 7 is the Kerr nonlinear parameter, which is posi- 
tive for our system (7 > 0). The nonlinear term can be 
discarded for low light intensities, for typical experimen- 
tal parameters see for example [ioj . 

Disorder is introduced to the lattice by randomly 
changing the width of each waveguide in a finite range 
W ± S where W is the mean value (typically 4/ira in our 
samples). The parameters j3 n and C n?n ±i can be calcu- 
lated numerically from the waveguides' width and from 
the separation between neighboring waveguides. As a re- 
sult of disorder the parameters j3 n become random in the 
range (3q ± A. By keeping the lattice periodic on average 
(the site's centers still have the lattice periodicity), the 
parameters C n?n ±i become independent of n to a very 
good approximation, meeting the conditions assumed by 
Anderson in his original model (diagonal disorder) [l|. A 
measure of disorder is then given by the ratio A/C ||. 

We now turn to the first set of experiments, designed 
to excite and measure individually all the localized eigen- 
modes of a disordered lattice, and to study the effect 
of nonlinear perturbations on them. For this purpose 
we have fabricated a lattice with N=99 sites and dis- 
order level A/C=l. The disordered lattice eigenmodes 
and eigenvalues can be calculated by diagonalizing a set 
of N equations (1) in the linear limit, substituting the 
calculated values of f3 n and C n?n ±i for this specific re- 
alization. Results are shown in Fig 2(a). The resulting 
band of eigenvalue deviates only slightly from the cosine- 
shaped band of a perfectly ordered lattice [19]. Several 
eigenmodes with eigenvalues near the edges of the band 



become exponentially localized in space, even though the 
density of states near the band edges is not significantly 
lower from the one at the band center [4(. Localized 
eigenmodes near the bottom of the band are flat phased, 
i.e. their wavefunction's amplitude is in-phase at all sites 
(see inset l)[2l|, while the localized eigenmodes at the 
top of the band are staggered, i.e. their wavefunction's 
amplitude has a tt phase flip between adjacent sites (inset 

2) . Notably, These localized eigenmodes are localized at 
well-separated regions in space, and have a similar width 
of about 10 lattice sites. Eigenmodes near the band cen- 
ter remain extended in the finite sized system (see inset 

3) . These modes will also be localized in an infinite sys- 
tem, but on a much larger length scale [22]. As disorder 
(A/C) is increased, a larger fraction of the modes be- 
comes localized within the finite lattice. 

In order to excite localized eigenmodes of this lattice 
we inject a wide beam (covering several lattice sites) at 
different locations across the lattice. At some input po- 
sitions, we observe significant expansion of the beam at 
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FIG. 3: (color online). The effect of weak nonlinearity on 
pure localized eigenmodes: cross-sections of the output light 
intensities (horizontal axis) at different power levels (verti- 
cal axis), showing that (a) flat phased localized modes tend 
to become more localized as nonlinearity is increased, while 
(b) staggered localized modes tend to delocalize. All cross- 
sections are normalized to unit maximum. 
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FIG. 4: (color online). The effect of disorder on wavepacket expansion in the linear and the nonlinear regimes. The results 
shown are normalized lattice averages of the output light distribution, when initially a single site is excited (see discussion in 
the text). A measure of localization is given by the mean square displacement from the input site a, and by the intensity at 
the origin site I . (a-d) Measurements in the linear case show the transition from ballistic transport to exponential (Anderson) 
localization as a function of disorder in: (a) A/C = 0, (b) A/C — 1, (c) A/C = 1.5, and (d) A/C = 3. Note the transition 
from ballistic transport in (a) to localization in (d) through the buildup of a central component and the suppression of the 
ballistic side lobes. Inset in (d) shows the localized distribution in semilog scale, demonstrating the exponential tails, (e)-(h): 
Measurements of the same lattices in the nonlinear case, showing that on average nonlinearity tends to increase localization for 
intermediate disorder levels (e.g. compare (b) and (f)). 



the output, similar in width to the expansion observed 
in a periodic lattice with the same average parameters 
(see fig 1(b) and 1(c)). At specific locations however, 
wave expansion is suppressed and prominent localization 
is evident (Fig 1(d)). In these cases, the input beam 
overlap significantly with a single localized eigenmode of 
the lattice, while the overlap with the other modes of the 
lattice is eliminated. Optimized localized output distri- 
bution were achieved using an input beam that covers 
about 10 lattice sites, corresponding to the typical width 
of the localized eigenmodes in this disorder level. 

While scanning the input position we observe several 
localized light distributions at the output. We compare 
the obtained localized intensity profiles to calculated lo- 
calized eigenmodes of the lattice in Fig 2(b). There is 
a clear correspondence between the experimentally ob- 
served localizations and the location and shape of all 
the calculated flat-phased localized eigenmodes associ- 
ated with the bottom of the eigenvalue band. To excite 
the staggered modes associated with the top of the band, 
the input beam was tilted with respect to the lattice to 
induce a 7r phase difference in the excitation of adjacent 
waveguides [l9| . The results of this excitation scheme are 
presented in Fig. 2(c). Again, a clear correspondence is 
found to the calculated localized staggered modes of the 
lattice. These results demonstrate the ability to excite 
pure Anderson localized eigenmodes. 

The effect of nonlinear perturbations on localized 
eigenmodes is studied by exciting a pure localized mode 
and increasing the input beam power. The intensities are 
kept way below those required to exhibit self-focusing in 



a periodic lattice with the same average parameters [20| , 
keeping the experiments in the weak nonlinear regime. 
Some localized modes are found to exhibit significant re- 
sponse to nonlinearity. The results of two such experi- 
ments are shown in Fig. 3, showing that weak positive 
nonlinearity tends to further localize flat-phased localized 
modes, but tends to de-localize staggered modes. 

These results can be understood on the basis of the the- 
ory developed in [12], which have shown that nonlinear 
shifts of a localized eigenmode's frequency (represented 
here by the parameter /3), can lead to derealization if 
the frequency crosses a resonance with other modes of 
the lattice. This condition can be satisfied in the case 
of weak disorder for the staggered localized modes, as 
nonlinearity shifts them from the edge of the band into 
the linear spectrum. Conversely, the flat-phased modes 
at the other edge of the band are shifted by nonlinearity 
out of the linear spectrum, thus they remain localized 
(see also discussion in [15[). 

We now turn to the second set of experiments, aimed 
to study the effect of disorder on wavepacket evolution in 
the linear and nonlinear regimes. This aspect is studied 
by injecting light into a single lattice site, thus exciting 
a tight (5-like wavepacket of all eigenmodes having non- 
vanishing overlap with the excited site. The wavepacket 
then evolves in the lattice, and the light distribution is 
measured at the output. We average the output patterns 
obtained, by separately exciting each site in the lattice 
while keeping the measurement-window centered about 
the input site position. The results of such measure- 
ments in the linear regime, taken in 5mm long samples 
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with different disorder levels, are shown in Fig. 4(a-d). 
Without disorder, single site excitation results in bal- 
listic propagation (wavepacket width grows linearly with 
time), recognized by a characteristic signature of two sep- 
arated lobes [1^, Ho] (Fig. 4(a)). At moderate disorder, 
a second component emerges, localized around the input 
site position (Fig. 4(b,c)). The localized and the ballis- 
tic components coexist in this regime. At high disorder 
a highly localized, exponentially decaying distribution is 
observed (Fig. 4(d)). This exponential decay of the ex- 
pansion profile is a hallmark of Anderson localization. 

These results offer a first direct look at the short time 
evolution of wavepacket s in ID disordered systems. It 
is known that for infinite disordered ID systems and for 
long time scales, wavepacket expansion is always fully 
suppressed. However, on short time scales, wavepackets 
do evolve 17, [23[ . The results in Fig. 4(a-d) reveal 
how localization is reached in the ID case; localization 
emerges from ballistic expansion through the continues 
buildup of a localized component and the suppression of 
a ballistic component. This dynamics is fundamen tally 
different from the one observed in 2D systems JjJ LL6[ , 
or that was indirectly measured in quasi- ID experiments 
(iq| . In these cases the expansion turns quickly from 
ballistic to diffusive, and becomes localized after much 
longer propagation times. In the ID case the diffusive 
dynamics is absent, as discussed for example in j23[. 

To study the effect of nonlinearity on wavepacket ex- 
pansion we repeat these measurements at increased light 
intensities. Again, we remain in the weak nonlinear 
regime. Results are shown in Fig. 4(e-h). On average, 
the results indicate increased localization at intermediate 
disorder levels; the intensity in the ballistic component 
is lower and the intensity in the localized component is 
higher under nonlinear conditions. This suggests that the 
buildup of the localized component and the suppression 
of the ballistic component happen faster under nonlinear 
conditions. This description holds for short time scales in 
which the ballistic component is still present. On much 
longer time scales, subdiffusive derealization due to non- 
linearity was predicted to take over [13[. 

In conclusion, we have directly studied localized eigen- 
modes and wavepacket expansion in disordered ID lat- 
tices that are described by the nonlinear version of the 
Anderson model. Two types of localized eigenmodes were 
measured, and nonlinear perturbations were shown to en- 
hances localization in one type and induce derealization 
in the other. The study of the expansion of wavepack- 
ets on short time scales in the presence of disorder has 
enabled a direct measurement of the transition from bal- 
listic wavepacket expansion to exponential localization. 
It was shown that in ID systems a ballistic and a lo- 
calized components can co-exist at intermediate times, 
while diffusive wavepacket expansion, observed in 2D and 
quasi- ID systems, is absent in the ID case. In the nonlin- 
ear regime, evidence is found for a faster transition into 



localization under nonlinear conditions. 
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